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Abstract. Free probability theory originated in the context of 

operator algebras, however, one of the main features of that theory 
is its connection with random matrices. Indeed, free probability 
can be considered as the theory providing concepts and notations, 
without relying on random matrices, for dealing with the limit 
N oo of N X iV-random matrices. 

One of the basic approaches to free probability, on which I will 
concentrate in this lecture, is of a combinatorial nature and centers 
axound so-called free cumulants. In the spirit of the above these 
arise as the combinatorics (in leading order) of N x A''-random 
matrices in the limit N = oo. These free cumulants are multi- 
linear functionals which are defined in combinatorial terms by a 
formula involving non-crossing partitions. 

I will present the basic definitions and properties of non-crossing 
partitions and free cumulants and outline its relations with freeness 
and random matrices. As examples, I will consider the problems 
of calculating the eigenvalue distribution of the sum of randomly 
rotated matrices and of the compression (upper left corner) of a 
randomly rotated matrix. 



1. Random matrices and freeness 

Free probability theory, due to Voiculcscu, originated in the context 
of operator algebras, however, one of the main features of that theory 
is its connection with random matrices. Indeed, free probability can 
be considered as the theory providing concepts and notations, without 
relying on random matrices, for dealing with the limit — > oo of 
A'" X iV-random matrices. 

Let us consider a sequence (v4jv)jveN of selfadjoint N x A^-random 
matrices A^. In which sense can we talk about the limit of these ma- 
trices? Of course, such a limit does not exist as a oo x oo-matrix and 
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there is no convergence in the usual topologies connected to opera- 
tors. What converges and survives in the limit are the moments of the 
random matrices. 

To talk about moments we need in addition to the random matrices 
also a state. This is given in a canonical way by the averaged trace: Let 
tr^v be the normahzed trace on N x A'"-matrices, i.e. for A — {aij)f^j^^ 
we have 

1 ^ 
1=1 

In the same way, we get the averaged trace tr^v E for N x A'^-random 
matrices, i.e. for A = {aij{uj))fj^i (where the entries random 
variables on some probability space fl equipped with a probability mea- 
sure P) we have 

1 ^ f 

tiN^HA) / au{uj)dP{uj). 

i=l "'^ 

Given these states trjv ® E, we can now talk about the k-th moment 
trjv ® HA%) of our random matrix A^, and it is well known that for 
nice random matrix ensembles these moments converge for TV — > oo. 
So let us denote by ak the limit of the A;-th moment, 

lim trAT ® E{A%) =: a*.. 

AT— »oo 

Thus we can say that the limit N = oo consists exactly of the collection 
of all these moments a^. But instead of talking about a collection of 
numbers we prefer to identify these numbers as moments of some 
variable A. Abstractly it is no problem to find such an A, we just take 
a free algebra A with generator A and define a state (p on Ahy setting 

cp{A'') := at. 

Of course, nothing deep has happened, this is just a shift in language, 
but it provides us with a more conceptual way of looking at the limit 
of our random matrices. Now we can say that our random matrices 
A]^ converge to the variable A in distribution (which just means that 
the moments oi An converge to the moments of A). We will denote 
this hy Aisf ^ A. Note that the nature of the limit A?" = oo is quite 
different from the case of finite N. In the latter case the An live in 
classical probability spaces of x A^-random matrices, whereas the 
N = oo limit object A is not of a classical nature any more, but lives 
in a 'non-classical probability space' given by some algebra A and a 
state (p. 
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1.1. Remark. One should note that for a selfadjoint operator A = A*, 
the collection of moments (or, equivalently, the state </? corresponding 
to these moments) corresponds also to a probability measure /ia on the 
real line, determined by 



(We can ignore the problem of non-uniqueness of this moment problem, 
because usually our operators A are bounded, which ensures unique- 
ness.) In particular, for a selfadjoint A^x A^-matrix A = A* this measure 
is given by the eigenvalue distribution of A, i.e. it puts mass 1/A^ on 
each of the eigenvalues of A (counted with multiphcity) : 



where Xi, . . . , Xn are the eigenvalues of A. In the same way, for a 
random matrix A, /ia is given by the averaged eigenvalue distribution 
of A. Thus, moments of random matrices with respect to the averaged 
trace tr^r ® E contain exactly that type of information in which one is 
usually interested when dealing with random matrices. 

1.2. Example. Let us consider the basic example of random matrix 
theory, expressed in our new language. Let Gn he the usual selfadjoint 
Gaussian N x A^-random matrices (i.e., entries above the diagonal are 
independently and normally distributed). Then the famous theorem of 
Wigner can be stated in our language in the form that 



where semi-circular just means that the measure is given by the 
semi-circular distribution (or, equivalently, the even moments of the 
even variable s are given by the Catalan numbers). 

Up to now, nothing crucial has happened, we have just shifted a bit 
the usual way of looking on things. A new and crucial concept, however, 
appears if we go over from the case of one variable to the case of more 
variables. Of course, again joint moments are the surviving quantities 
in multi-matrix models (even if it is now not so clear any more how to 
prove this convergence in concrete models) and we can adapt our way 
of looking on things to this situation by making the following definition. 





i=l 



where s is a semi-circular variable. 
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1.3. Definition. Consider N x A^-random matrices Aj^^ , . . . , and 
variables Ai, . . . , Am (living in some abstract algebra A equipped with 
a state </?) . We say that 

{A^^\ • • • ' ^jT^) {Ai,---, Am) in distribution, 

if 

hm tr^ ® E[^(;^) • • • A^^''] = (^(Ai • • • AJ 

N—too 

for all choices oi k, 1 < ii, . . . ,ik < m. 

1.4. Remark. The Ai, . . . , A^, arising in the limit of random matrices 
are a priori abstract elements in some algebra A, but it is good to know 
that in many cases they can also be concretely realized by some kind 
of creation and annihilation operators on a full Fock space. Indeed, 
free probability theory was introduced by Voiculescu for investigating 
the structure of special operator algebras generated by these type of 
operators. In the beginning, free probability had nothing to do with 
random matrices. 



1.5. Example. Let us now consider the example of two independent 
Gaussian random matrices G^^\ G^^ (i.e., each of them is a selfadjoint 
Gaussian random matrix and all entries of G^^^ are independent from 
all entries of G^^). Then one knows that all joint moments converge, 

and we can say that (G^^G^^^) — >• (si,S2), where Wigner tells us 
that both Si and S2 are semi-circular. The question is: What is the 
relation between Si and S2? Does the independence between G^^"* and 
G^^ survive in some form also in the limit? The answer is yes and is 
provided by a basic theorem of Voiculescu which says that si and S2 
are free in the following sense. 

1.6. Definition. Let ^ be a unital algebra and (p : A ^ C a, linear 
functional on A, which is unital, i.e., ip{l) — 1. Then oi, . . . , G A 
are called free (with respect to </?) if 

<)^bi(«i(i)) • • •Pfc(«i(fe))] = 

whenever 

• pi, . . . ,pk arc polynomials in one variable 

• i{2) ^ i(3) • • • i{k) (only neigbouring elements are 
required to be distinct) 

• (fi\pj{ai(^j)\ = for all j = 1, . . . , /c 
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1.7. Remark. Note that the definition of freeness can be considered 
as a way of organizing the information about all joint moments of 
free variables in a systematic and conceptual way. Indeed, the above 
definition allows to calculate mixed moments of free variables in terms 
of moments of the single variables. For example, if a, b are free, then 
the definiton of freeness requires that 

(^[(a-(^(a)-l)(6-^(6)-l)] = 0, 

which implies that 

ip{ab) = (p{a) ■ ip{b) if a, b are free. 
In the same way, 

cp[{a - <p{a) ■ l){b - <p{b) ■ l)(a - <p{a) •!)(&- cp{b) • 1)] = 
leads finally to 

ip{abab) = ip{aa) ■ ip{b) ■ ip{b) + ip{a) ■(p{a) ■ (p{bb) —(p{a) -ifib) ■ ip{a) ■ (p{b). 

Analogously, all mixed moments can (at least in principle) be calculated 
by reducing them to alternating products of centered variables as in the 
definition of freeness. 

Thus the statement 'si, S2 are free and each of them is semicircular' 
determines all joint moments in si and S2- 

Formulating our knowledge about the joint moments of Si and S2 in 
this peculiar way might look not very illuminating on first sight, but it 
will turn out that recognizing this notion of freeness as the organizing 
principle for the collection of moments adds a new perspective on the 
limit of random matrices. 

In particular, we are now in the context of non-commutative prob- 
ability theory which consists mainly of the doctrine that one should 
use notations and ideas from classical probability theory in order to 
understand problems about non- commutative algebras. 

Free probability theory can be described as that part of non-com- 
mutative probability theory where the notion of 'freeness' plays an es- 
sential role. Furthermore, according to the basic philosophy of Voiculescu 
this notion of freeness should be considered (and investigated) in anal- 
ogy with the classical notion of 'independence' - both freeness and in- 
dependence prescribe special relations between joint moments of some 
variables. (Of course, both cases correspond to very special, but also 
very fundamental situations.) 
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One of the most interesting features of freeness is that this concept 
appears in at least two totally different mathematical situations. Orig- 
inally it was introduced by Voiculescu in the context of operator alge- 
bras, later it turned out that there is also some relation, as described 
above, with random matrices. This gives a non-trivial connection be- 
tween these two different fields. For example, modelling operator alge- 
bras by random matrices has led to some of the most impressive results 
about operator algebras in the last years. 

Furthermore, apart from the concrete manifestations of freeness via 
random matrices or operator algebras, there exist also an abstract prob- 
abilistic theory of freeness, which shows the depth of this concept and 
which I want to address in the following. 



2. Combinatorial approach to free probability: 
non-crossing partitions and free cumulants 

'Freeness' of random variables is defined in terms of mixed moments; 
namely the defining property is that very special moments (alternating 
and centered ones) have to vanish. This requirement is not easy to 
handle in concrete calculations. Thus we will present here another ap- 
proach to freeness, more combinatorial in nature, which puts the main 
emphasis on so called 'free cumulants'. These are some polynomials in 
the moments which behave much better with respect to freeness than 
the moments. The nomenclature comes from classical probability the- 
ory where corresponding objects are also well known and are usually 
called 'cumulants' or 'semi-invariants'. There exists a combinatorial 
description of these classical cumulants, which depends on partitions 
of sets. In the same way, free cumulants can also be described com- 
binatorially, the only difference to the classical case is that one has to 
replace all partitions by so called 'non-crossing partitions'. 

This combinatorial description of freeness is due to me [8, 9] (see 
also [3]); in a series of joint papers with A. Nica [4, 5, 6] it was pursued 
very far and yielded a lot of new results in free probability theory. For 
more information on other aspects of freeness, in particular the original 
analytical approach of Voiculescu, one should consult the papers [10, 
11, 13], the collection of various articles [13], or the monographs [14, 1] 

2.1. Definitions. A partition of the set S := {1, . . . , n} is a decompo- 
sition 



n^{V^,...,Vr} 
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of S into disjoint and non-empty sets Vi, i.e. 

r 

VinVj^$ (i,j^l,...,r;i^j) and S ^[jVi. 



i=l 



We call the Vi the blocks of tt. 
For 1 < p,q < n we write 

p '-^■K q if p and q belong to the same block of tt. 

A partition tt is called non-crossing if the following does not occur: 
There exist 1 < Pi < < P2 < ?2 < ^ with 

The set of all non-crossing partitions of {1, .... n} is denoted by NC(n). 

Non-crossing partitions were introduced by Kreweras [2] in a purely 
combinatorial context without any reference to probability theory. 

2.2. Examples. We will also use a graphical notation for our parti- 
tions; the term 'non-crossing' will become evident in such a notation. 
Let 

5 = {1,2,3,4,5}. 

Then the partition 



7r = {(l,3,5),(2),(4)} 
is non-crossing, whereas 



2345 
I I I I 



7r = {(l,3,5),(2,4)} 



2345 



is crossmg. 



2.3. Remcirks. 1) In an analogous way, non-crossing partitions NC{S) 
can be defined for any linearly ordered set S; of course, we have 

NC{Si) ^ NC{S2) if = ^82. 

2) In most cases the following recursive description of non-crossing 
partitions is of great use: a partition tt ist non-crossing if and only if 
at least one block K G tt is an interval and Tr\V is non-crossing; i.e. 

V E 71 has the form 

V — {k, k + 1, . . . ,k + p) for some 1 < k <n and p>0,k-\-p<n 
and we have 

7r\V e NC(1, ...,k-l,k+p + l,...,n)^ NC(n -(p+ 1)). 
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Example: The partition 



23456789 
I U I I U I 







{(1,10), (2, 5, 9), (3, 4), (6), (7, 8)} = 
can, by successive removal of intervals, be reduced to 
{(1,10), (2,5,9)}={(1,5), (2,3,4)} 

and finally to 

{(1,5)}={(1,2)}. 

3) By writing a partition vr in the form tt = {Vi, . . . , Vr} we will always 
assume that the elements within each block Vi are ordered in increasing 
order. 

2.4. Definition. Let {A, 9?) be a probability space, i.e. ^ is a unital 
algebra and 9? : ^ — > C is a unital linear functional. We define the 
(free or non-crossing) cumulants 

kn-.A^'^C (neN) 

(indirectly) by the following system of equations: 

(fi{ai...an)= ^ K[ai,...,an\ (oi, . . . , a„ e ^), 

neNC{n) 

where denotes a product of cumulants according to the block struc- 
ture of tt: 

KWi, . . . , a„] := /cyjai, . . . , a„] . . . kv,[ai, . . . , a„] 

ior n ^ {Vi, . . . ,Vr} e NC{n) 

and 

/i;y[ai,...,a„] := A;#y(a^i,...,a„J for V ^ {vi, . . . ,vi). 

2.5. Remarks and Examples. 1) Note: the above equations have 
the form 

ip{ai . . . ttn) = kn{ai, . . . , a„) + smaller order terms 

and thus they can be resolved for the kn{ai, . . . , an) in a unique way. 
2) Examples: 

• n — 1 

(/?(ai) = k\[ai] = A;i(ai), 

thus 

ki{ai) = </?(ai). 
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• n = 2 

<y5(«i«2) = ku[ai, 0.2] + k\ i[ai, 02] 

= k2{ai, 02) + ki{ai)ki{a2), 

thus 

/C2(ai,a2) = (p{aia2) - 99(01)99(02). 

• n = 3 

9?(ai02a3) = k\_u[ai, (12. as] + /ci u [«■!! «'2, o-s] + i[ai,«2,a3] 
+ A;uj [01,02,03] + /ci I 1(01,02,03] 
= A;3(oi, 02, 03) + /ci (01)^2(02, 03) + k2{ai, a2)ki{a3) 
+ k2{ai,as)ki{a2) + /ki(oi)/i;i (02)^:1(03), 

and thus 

A;3(ai, 02, 03) = 93(010203) - 99(01)99(0203) - 99(0103)9^3(02) 
- 99(0102)99(03) + 299(01)99(02)99(03). 

3) For n = 4 we consider the special case where all <^(oj) = 0. Then 
we have 

^4(01, 02, 03, 04) = 9^5(01020304) - 99(0102)99(0304) - 99(0104)92(0203). 

4) The kn are multi-linear functionals in their n arguments. 

The meaning of the concept 'cumulants' for freeness is shown by the 
following theorem. 

2.6. Theorem. Consider oi,...,Oto G A. Then the following two 
statements are equivalent: 

i) oi, . . . , ttm are free. 

ii) mixed cumulants vanish, i.e.: We have for all n > 2 and for all 
1 < i(l), . . . , i{n) < m: 

kn{ai(i), ai{n)) = 0, 
whenever there exist 1 < p,q < n with i(p) i{q)- 

2.7. Remcirks. 1) An example of the vanishing of mixed cumulants is 
that for a,b free we have k3{a,a,b) — 0, which, by the definition of ks 
just means that 

(p{aab) — (p{a)(p{ab) — (p{aa)(p{b) — (p{ab)(p{a) + 2(p{a)ip{a)ip{b) — 0. 

This vanishing of mixed cumulants in free variables is of course just a 
reorganization of the information about joint moments of free variables 
- but in a form which is much more useful for many applications. 
2) The above characterization of freeness in terms of cumulants is the 
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translation of the definition of freeness in terms of moments - by using 
the relation between moments and cumulants from Definition 2.4. One 
should note that in contrast to the characterization in terms of moments 
we do not require that i(l) 7^ i{2) 7^ ■ ■ ■ 7^ i{n) or ip{ai) = 0. (That's 
exactly the main part of the proof of that theorem: to show that on 
the level of cumulants the assumption 'centered' is not needed and that 
'alternating' can be weakened to 'mixed'.) Hence the characterization 
of freeness in terms of cumulants is much easier to use in concrete 
calculations. 



3. Addition of free veiriables 
3.1. Notation. For a random variable a e ^ we put 

and call {k^)n>i the (free) cumulants of a. 

Our main theorem on the vanishing of mixed cumulants in free vari- 
ables specialises in this one-dimensional case to the linearity of the 
cumulants. 



3.2. Proposition. Let a and b be free. Then we have 
k^+f'^k^ + k^^ foralln>l. 

Proof. We have 

k^+'' = kn{a + b,...,a + b) 

= kn{a, ...,a) + knih, . . . , 6) 
— h"- j^h^ 



because cumulants which have both a and h as arguments vanish by 
Theorem 2.6. □ 



Thus, the addition of free random variables is easy to describe on the 
level of cumulants; the cumulants are additive under this operation. It 
remains to make the connection between moments and cumulants as 
explicit as possible. On a combinatorial level, our definition specializes 
in the one-dimensional case to the following relation. 
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3.3. Proposition. Let (m„)n>i and {kn)n>i be the moments and free 
cumulants, respectively, of some random variable. The connection be- 
tween these two sequences of numbers is given by 

neNC{n) 

where 

:= k^Vi ■ ■ ■ k#Vr for tt = {T^, . . . , K}- 
Example: For n = 3 we have 

ms = kuj +k\ij+ku\+ ^lli + i i 
= /ca + 3kik2 + kl- 

For concrete calculations, however, one would prefer to have a more 
analytical description of the relation between moments and cumulants. 
This can be achieved by translating the above relation to corresponding 
formal power series. 

3.4. Theorem. Let {mn)n>i and (/c„)n>i be two sequences of complex 
numbers and consider the corresponding formal power series 



M{z) := l + ^m^^^^, 

n=l 

oo 



n=l 

Then the following three statements are equivalent: 

(i) We have for all n e N 

rrin = E /CjT = k#Vi ■ ■ ■ k#Vr 

neNC{n) n={Vi,...,Vr}eNC{n) 

(ii) We have for all n e N (where we put mo := 1) 



(iii) We have 



S=l <i,...,ise{0,l,...,n— s} 
2lH +is=n—s 



C[zM{z)] = M{z). 
Proof. We rewrite the sum 
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in the way that we fix the first block Vi of tt (i.e. that block which 
contains the element 1) and sum over all possibilities for the other 
blocks; in the end we sum over Vi: 



n 

rrin 



S=l Vi with #Vi = S TreNC(n) 

where tt = { Vi , . . . } 

If 

Vi = {Vi = 1,V2, . . .,Vs), 

then TT = {Vi, . . . } e NC{n) can only connect elements lying between 
some Vk and Vk+i, i.e. vr = {Vi, V2, . . . , K} such that we have for all 
J = 2, . . . , r: there exists a k with Vk <Vj < t'fe+i. There we put 

fs+i := n + 1. 

Hence such a tt decomposes as 

TT = U TTi U • • • U 7fs, 

where 

iTj is a non-crossing partition of {vj + 1, + 2, . . . , Vj+i — 1}. 
For such TT we have 

/Cjr — k^y^kj^i ■ ■ ■ kj^^ — kgkj^^ . . . kj^^, 
and thus we obtain 



n 



S=l l=V\<V2<---<Vs<n •n-=ViUiriU---U7rs 

frjeJVC(v3+l,-.-,fj + l-l) 

n 

= E*. E ( E E 

s=l l=?;i<?;2<---<i)s<T» 7rieA''C(?;i+l,...,t)2-l) ifs6ArC(?;3+l,...,n) 

n 

— ^ ^ ^ ^ '^V2—v\ — l ■ ■ ■ T^n—Vs 

s=l \=V'i<V2<- ■ ■ <v s<n 
n 

S=l ii,...,ise{0,l,...,n-3} 
ilH -i-is+s=n 

This yields the implication (i) =^ (ii). 

We can now rewrite (ii) in terms of the corresponding formal power 
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series in the following way (where we put mo := ko := 1): 

oo 

M(z) = l + ^z'^m^ 



n=l 

oo n 



n=0 



n=l S=l ii,...,ise{0,l,...,n— s} 

«lH — +is=n—s 

OO oo 
s=l i=0 

= C[^M(^)]. 
This yields (iii). 

Since (iii) describes uniquely a fixed relation between the numbers 
{kn)n>i and the numbers (m„)„>i, this has to be the relation (i). □ 

If we rewrite the above relation between the formal power series in 
terms of the Cauchy transform 

oo 

G{z) := Y: 
and the i?-transform 

oo 

R{z) ■.= Y,kn+iz\ 

n=0 

then we obtain Voiculescu's formula. 

3.5. Corolleiry. The relation between the Cauchy transform G{z) and 
the i?-transform R{z) of a random variable is given by 

G[R{z) + \]=z. 

Proof. We just have to note that the formal power series M{z) and 
C{z) from Theorem 3.4 and G(2;), i?(2;), and K{z) = R{z) + ^ are 
related by: 

G{z) = -Mil) 
z z 

and 

C{z) = 1 + zR{z) = zK{z), thus K{z) = 

z 

This gives 

A-1G(.)] ^ ^C[G(.-)1 ^ g^C[iM(l)] ^ ^M(i) ^ 
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thus K[G{z)] = z and hence also 

G[R{z) + -] = G[K{z)\ = z. 
z 

□ 

3.6. Free convolution. The above results give us a quite effective 
tool for calculating the distribution of the sum a + 6 of free variables 
from the distribution of a and the distribution of h. In analogy with 

the usual convolution (which corresponds to the sum of independent 
random variables) wc introduce the notion ffl of /ree convolution as 
operation on probability measures by 

fia+b = A*a ffl A'-b if ^L, b are free. 

Then we know that free cumulants and the i?-transform linearize this 
free convolution. 

In particular, we also have the free convolution powers 

/X™'" := /X ffl • • • ffl /X (r-times) 

of n, which are in terms of cumulants characterized by 

If we are given free variables a and b and we want to calculate the 
distribution of a + 6, then we calculate the i?-transforms Ra and Rb 
and get thus by the linearization property the i?-transform of a + 6, 

Ra+b = Ra + Rb- 

It remains to extract the distribution out of this. From the /?-transform 
we can get the Cauchy transform Ga+b by Corollary 3.5, and then we 
use the classical Stietljes inversion formula for extracting the distri- 
bution from this. In general, the relation between i?-transform and 
Cauchy transform might lead to equations which have no analytic so- 
lutions, however, in many concrete case these equations can be solved. 
For example, if we put = |(5o + ^i), then the above machinery shows 
that the distribution ffl is given by the arcsine law. 

3.7. Application: random sum of matrices [7]. Fix two (com- 
pactly supported) probability measures n and v on the real line and 
consider deterministic (e.g., diagonal) N x iV-matrices and Cn, 
whose eigenvalue distribution converges, for ^ oo, towards the given 
measures. To put it in the language introduced in Section 1, we assume 
that 

Af^ — > a with /la — IJ' 
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and 

Cjv — > c with /ic — T^- 
Now we rotate Cjv against An randomly by replacing Cjv by 

Bn '■— UnCnUn: 

where C/at is a random Haar unitary matrix from the ensemble of uni- 
tary X iV-matrices equipped with the Haar measure. Of course, the 
eigenvalue distribution of Bjsj is the same as the one of C]v, however, 
any definite relation between the cigcnspaccs of Ajy and the eigenspaces 
of Cjv has now been destroyed. Aj^ and Bp^ are in the limit N ^ oo 
generic realizations of the given eigenvalue distributions /i and i/. The 
question which we want to address is: What is the eigenvalue distribu- 
tion of the sum An + Bn in the limit N —>■ oo, i.e. what can we say 
about 

An + Bn ^? 

A version of the theorem of Voiculescu about the connection between 
random matrices and freeness tells us that An and Bn become free in 
the limit N oo, i.e. it yields that 

{An - Bn) — > {a, b) with /Xq = //, fih = z^, and a, b free. 

Thus we know that the eigenvalue distribution of An + Bn converges 
towards the distribution of a + 6 where a and b are free. But the distri- 
bution of a -|- 6 can be calculated with our tools from free probability 
in a very effective and systematic way by using the i?-transform ma- 
chinery. For example, if we take the generic sum of two projections of 
trace 1/2, (i.e., ji = v = |((5o + 5i), then our example from above shows 
us that the distribution of 



/I 



/I 



N 



u 



\ 



\ 



■I 



is, in the limit N ^ oo, given by the arcsine law. 

4. Multiplication of free variables 

Finally, to show that our description of freeness in terms of cumulants 
has also a significance apart from dealing with additive free convolution, 
we will apply it to the problem of the product of free random variables: 
Consider oi, . . . , a„, 6i, . . . , 6„ such that {oi, . . . , a„} and {bi, . . . , bn} 
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are free. We want to express the distribution of the random variables 
oi&i, . . . , a„6„ in terms of the distribution of the a's and of the 6's. 

4.1. Notation. 1) Analogously to A;^ we define for 



where 

(/7v^[ai,...,a„] := (/7(a„i ...a„J for V ^ {vi, . . . ,vi). 
Examples: 



</?i_u [01,02,03] = (p{aia2a3) 
<P\ u [oi, a2, 03] = V'(ai)v?(a2a3) 
I [oi, 02, 03] = V9(aia2)v^(a3) 
[01,02,03] = V9(aia3)v9(a2) 

(/?| I |[oi, O2, O3] = (^(0i)(^(02)(^(03) 



2) Let (7, TT e NC{n). Then we write 



if each block of a is contained as a whole in some block of tt, i.e. a can 
be obtained out of tt by refinement of the block structure. 
Example: 



{(1),(2,4),(3),(5,6)}<{(1,5,6),(2,3,4)} 



With these notations we can generalize the relation 



n = {V^,...,Vr}eNC{n) 



the expression 



, an] Vvi[ai, . . . , o„] . . . ipvAai, a^], 



a < TT 





TTeNC{n) 



in the following way. 





TveNC(n) 
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Consider now 

{ai, . . . ,a„}, {61, . . . free. 

We want to express alternating moments in a and b in terms of moments 
of a and moments of b. We have 

ip{aibia2b2 . . . anbn) = ^ K[ai,bi,a2,b2, . . . ,an,bn\. 

n£NC(2n) 

Since the a's are free from the 6's, Theorem 2.6 tells us that only such 
TT contribute to the sum whose blocks do not connect a's with 6's. But 
this means that such a tt has to decompose as 

7r = 7riU7r2 where tti e A'"C(1, 3, 5, . . . , 2n — 1) 

TTs e NC{2,A,6,...,2n). 

Thus we have 

(^(01610262 . . . anbn) 

— ^7ri[ai) ^2, • • • , On] ■ ^7r2[^l) ^2, • • • , ^n] 

TTj g]VC(odd),ir2 6iVC(euen) 
7r^U7r2eiVC(2n) 

= ^ ^A;^Jai,a2,...,a„] • ^ A;^2[^i> ^2, ■ ■ ■ , &n]) ■ 

Note now that for a fixed tti there exists a maximal element a with 
the property tti U cr e NC {2n) and that the second sum is running over 
all 7r2 < (7. 

4.2. Definition. Let tt G NC{n) be a non-crossing partition of the 
numbers 1, . . . ,n. Introduce additional numbers 1, . . . ,n, with alter- 
nating order between the old and the new ones, i.e. we order them in 
the way 

1122... nn. 

We define the complement K{t:) of tt as the maximal a e NC{1, . . . ,n) 
with the property 

nUa e NC(l,l,...,n,n). 

If we present the partition tt graphically by connecting the blocks in 
1, . . . ,n, then a is given by connecting as much as possible the numbers 
1, . . . ,n without getting crossings among themselves and with tt. 

(This natural notation of the complement of a non-crossing partition 
is also due to Kreweras [2]. Note that there is no analogue of this for 
the case of all partitions.) 



18 



ROLAND SPEICHER 



With this definition we can continue our above calculation as follows: 
ip{aibia2b2 ■ ■ ■ anbn) 

= X] (^7ri[ai,a2,---,an] • KjpiM, ■ ■ ■ ,bn]^ 



■KieNC(n) 



= XI ^7ri[ai,a2,---,an] •</?^^(7rl)[^'l,^'2,■■■,&n]■ 
7^l6^^C(r^) 

Thus we have proved the following result. 

4.3. Theorem. Consider 

{ai, . . . ,a„}, {61, . . . free. 

Then we have 

(/7(ai6ia2&2 ■ ■ ■ OnM = XI ^'^t 

neNC{n) 

Similar to the additive case one can translate this combinatorial de- 
scription of the product of free variables in an analytic way in terms 
of the so-called ^-transfrom. However, this is more complicated as in 
the case of the i?-transform and we will not address this problem here. 
Instead, we want to show that the above combinatorial description of 
the product of free variables can lead to quite explicit (and unexpected) 
results without running through an analytic reformulation. Such a re- 
sult is given in our final application to the problem of the compression 
of a random matrix. 

4.4. Application: Compression of random matrix. Consider, as 
in Section 3.7, a sequence of deterministic x A^-matrices Cat with 
prescribed eigenvalue distribution /i in the limit — > 00 and consider 

the randomly rotated version Ajy := UnCnU^]- of this matrix. The 
question we want to address is the following: Can we calculate the 
eigenvalue distribution of upper left corners of the matrices A^. For- 
mally, we get these corners by compressing Ajv with projections of the 
form 

/I \ 



■'AT •-- 
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where tr7v(-P/v) — ^ a for some fixed a with < a < 1. Thus we 
ask for the eigenvalue distribution of PnA^Pn in the hmit N ^ oo. 
(However, we have to calculate this in the compressed space, throwing 
away the bunch of trivial zero eigenvalues outside of the non-trivial 
corner of PatAtv-P/v, i.e., we consider P^AnPn not as x iV-matrix, 
but as aN x ctA^'-matrix.) 

Now, again by Voiculescu's theorem about asymptotic freeness of 
random matrices, we know that 

(^iv,-Pjv) (a,p), 

where a has the prescribed distribution p is a projection of trace 
a and a and p are free. Thus, the answer for our question on the 
distribution of corners in randomly rotated matrices is provided by 
calculating the distribution of pap in the compressed space, i.e. by 
calculating the renormalized moments 

a 

which is, by the trace property of (f and the projection property = p, 
the same as 

yiiapn 

a 

This fits now exactly in the above frame of calculating the moments of 
products of free variables, in the special case where the second variable 
is a projection of trace a. Using p^ — p for all A; > 1 and ip{p) — a 
gives 

Vk{^)[p,P, = (p{p...p)^{p...p)--- = a'-^^^^l, 

where |-f^(7r)| denotes the number of blocks of K(n). Wc can express 
this number of blocks also in terms of n, since we always have the 
relation 

|7r| + |X(7r)| + 
Thus we can continue our calculation of Theorem 4.3 in this case as 

-^[(ap)"] = - V A;.[a,...,a]a"+H-I 
a a ^-^ 

■K&NC(n) 

" X] -j;^^k^[aa,...,aa], 

n£NC(n) 

which shows that 

knipap, . . . ,pap) — —kn(cKa, . . . , aa) 
a 
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for all n. By our remarks on the additive free convolution, this gives 
the surprising result that the renormalized distribution of pap is given 

by 

In particular, for a = 1/2, we have 

ffl2 m 
^J'pap — /^l/2a — /^l/2o m /^l/2a- 

This means that the distribution of the upper left corner of size 1/2 of 
a randomly rotated matrix is, apart from rescaling with the factor 1/2, 

the same as the distribution of the sum of the considered matrix and 
another randomly rotated copy of itself. E.g., if wc take the example 
^ — + 5i), then the corner of size 1/2 of such a randomly rotated 
projection has as eigenvalue distribution the arcsine law. 
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